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Abstract 

Harmonic and anharmonic vibrational frequencies, as well as infrared and Raman 
intensities, are calculated for MgS0 4 nH 2 0 (n=l-3). Electronic structure theory at the second 

order Mpller-Plesset perturbation theory (MP2) level with a triple-zeta + polarization (TZP) 
basis set is used to determine the geometry, properties, and vibrational spectra of pure and 
hydrated MgSC>4 salts. The direct vibrational self-consistent field (VSCF) method and its 
correlation corrected (CC-VSCF) extension are used to determine anharmonic corrections to 
vibrational frequencies and intensities for the pure MgS04 and its complex with one water 
molecule. Very significant differences are found between vibrational spectra of water 
molecules in complexes with MgS0 4 and pure water. Some of the O-H stretching 

frequencies are shifted to the red very significantly (by up to 1500-2000 cm ') upon 
complexation with magnesium sulfate. They should be observed between 1700 and 3000 
cm ', in a region very different from the corresponding O-H stretch frequency region of pure 

water (3700-3800 cm" 1 ). In addition, the S0 2 stretching vibrations are found at lower 

frequency regions than the water vibrations. They can serve as unique identifiers for the 
presence of sulfate salts. The predicted infrared and Raman spectra should be of valuable 
help in the design of future missions and analysis of observed data from the ice surface of 
Jupiter’s moon Europa that possibly contains hydrated MgSQ 4 salts. 
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Introduction 

Recent observations from the Galileo satellite indicate that three of the Jupiter 

1-4 

moons, Europa, Ganymede, and Callisto, may have subsurface oceans . Because water 
is an essential ingredient of life, possible existence of the subsurface ocean and the nature 
of its composition are of great interest to planetary scientists and astrobiologists. Of the 
three moons, Europa has attracted the most attention. Data from Galileo’s near-infrared 
mapping spectrometer (NIMS) indicate the possibility of hydrated salts on Europa’s 
surface ' . The data suggested that Europa’s icy crust might contain hydrated sulfates of 
magnesium and sodium. However, the analysis of NIMS spectra is mainly based on the 
small shifts and broadening of the water absorption spectra in hydrated salts. Thus, the 
interpretation of NIMS spectra is not straightforward and does not provide a unique 

7 

assignment of the observed bands . On the other hand, a definitive identification of the 
existence of hydrated salts on Europa is of great interest to astrobiologists who have ranked 
Europa only behind Mars as a potential candidate for extraterrestrial life. Laboratory studies 
of hydrated MgS0 4 and Na 2 S0 4 salt minerals under low-temperature and ultrahigh vacuum 
conditions similar to those present on the surface of Europa showed that such salts could 

g 

exist on the surface of Jupiter’s satellite for geologic timescales . 

In this study, we use ab initio quantum chemistry calculations to obtain accurate 
vibrational spectra of magnesium sulfate, MgS0 4 , hydrated with one, two, and three water 
molecules. Our calculations provide positions and intensities of vibrational bands for both 
MgS0 4 and water molecules in the hydrated complexes that can serve as unique identifiers 
for the presence of such salts. The results of our theoretical calculations can assist in the 
analysis and interpretation of the experimentally observed spectra, and in the design of 
better detectors for future observations of the Jupiter moon surfaces. 

Theoretical Approach 

Stationary points on the potential energy surfaces of magnesium sulfate and its 
complexes with one, two, and three water molecules are located using second-order 
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M0ller-Plesset perturbation theory (MP2) 9 with a triple-zeta + polarization (TZP) basis set: 
Dunning's basis set for O and H atoms 10 and the McLean-Chandler basis set for Mg and 
S 11 . This level of theory is also used to calculate harmonic and anharmonic vibrational 
frequencies. Stationary points (minima and transition states) are determined by optimization 
using analytic gradients of MP2/TZP energies. Second derivative (Hessian) matrices are 
calculated numerically using double differencing of analytical gradients. Diagonalization of 
the Hessian matrices provides harmonic frequencies and normal mode coordinates. 
Anharmonic vibrational frequencies are computed for the smaller systems: pure MgS0 4 and 
its complex with one water molecule. Anharmonic corrections are obtained using the 
vibrational self-consistent field (VSCF) method 12 ’ 13 and its correlation corrected (CC- 
VSCF) extension via second-order perturbation theory 14 . The calculations are performed 
using the combined ab initio / CC-VSCF approach, where anharmonic vibrational 
frequencies are computed directly from an ab initio electronic structure program. Such 
direct techniques that employ ab initio potential energy surfaces have been developed 
recently 15 ’ 16 and have very important advantages over the techniques that use analytical 
potential functions. They allow the application of anharmonic vibrational spectroscopy 
methods to any molecular system of moderate size (up to 10-15 atoms), for which 
analytical potential functions are not available 15 . The VSCF method is used to solve the 
vibrational Schrodinger equation to obtain vibrational wavefunctions and energies. In the 
VSCF approximation, the total vibrational wavefunction in mass-weighted normal 
coordinates is represented as a product of one-dimensional functions (each depending on 
one normal coordinate), and the resulting one-dimensional VSCF equations are solved self- 
consistently. The VSCF energies are further corrected for correlation effects between the 

vibrational modes using second order perturbation theory (CC-VSCF) 14 . To make the 
integrals involved in the CC-VSCF calculations more tractable and to reduce the number of 
electronic structure computations required, a pairwise coupling approximation is used for 


4 


the potential in the normal-mode representation 14 . In this approximation, the potential of 
the system is represented by the sum of separable (single mode) terms and pair coupling 
terms, neglecting triple coupling of normal modes and higher-order interactions. This 
approximation was shown to work well in a number of previous calculations of hydrogen- 

bonded systems 17 . A grid representation for the potential energy terms is used, where one- 
dimensional separable terms for each normal mode are represented by a grid of 16 points, 
while pair coupling potential terms for each pair of normal modes is represented by a 

16x16 square grid. The potential energy values are computed on less dense (8x8) grids 

directly from the ab initio program and subsequently interpolated into 16x16 grids using 

polynomial interpolation. The number of potential energy points that needs to be computed 

a/ *(n _n » 

in the above algorithm is equal to A points = iV mode * JV grid + — 1 — — * N gnd where 

Nmode is the number of normal modes of the system; N gr id is the number of grid points in 
one dimension. The same MP2/TZP ab initio level used to obtain equilibrium structures and 
harmonic normal modes is also used to compute potential energies on grids. Full details of 
the CC-VSCF method combined with an ab initio electronic structure code is given in Ref. 
[15, 17]. All calculations in this study are performed using the electronic structure package 

GAMESS 18 . 

Previous studies showed that the MP2/TZP level of theory gives satisfactory results 
for geometrical structures and spectroscopic properties of different hydrogen-bonded 
complexes: water clusters 17 , complexes of negative and positive ions with water 17 , and 
complexes of acids (HC1 and HF) with water 19 ’ 2 ®. The combined ab initio / CC-VSCF 
technique that accounts for anharmonicities and couplings between vibrational modes (in 
conjunction with the MP2/TZP ab initio potentials) was shown to predict anharmonic 
frequencies for the highest frequency stretching vibrations with the accuracy of 30-50 cm 1 
compared with available experimental data 17 * 19 ' 20 . Though this is not at the level of 
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spectroscopic accuracy, it is a substantial improvement over the harmonic level theoretical 
predictions. Furthermore, even at this level of accuracy theoretical calculations can be used 
to help analyze observed spectra and predict important spectroscopic features. The method 
has been shown to have sufficient accuracy even for highly anharmonic hydrogen-bonded 
systems. Complexes of magnesium sulfate with water considered in this study are expected 
to be very anharmonic. Reliable prediction of their spectroscopic properties requires 
methods that go beyond the harmonic approximation. The combined CC-VSCF / ab initio 
method should give useful predictions for the most important spectroscopic features 
characteristic for such hydrated magnesium sulfate systems. 

Results and Discussion 

A. Equilibrium geometries, energetics, and charge distribution 

The structure of the pure MgS0 4 and pure water molecule obtained at the MP2/TZP 
level of theory are shown in Fig. 1 . Also shown are partial charges, obtained using the 

natural bond orbital (NBO) analysis 21 . The figure shows large positive charges on S and 
Mg atoms. The smallest 1:1 complex (MgS0 4 H 2 0) is found to have two stable isomers, 

shown in Fig. 2. The one that corresponds to the global minimum has a hydrogen bond 
between the hydrogen atom of the water molecule and one of the oxygens of \lgSO4- The 
O-H bond that participates in this hydrogen bonding is significantly elongated compared to 
O-H distances of the pure water molecule (by 0.086 A). One of the S=0 bonds is also 
elongated by 0.056 A, while the other S=0 bond is slightly shortened (by 0.01 A). In 
addition to this very strong hydrogen bond, there is also a strong electrostatic interaction 
between Mg atom and the oxygen of H2O, also present in the other isomer (Fig. 2b). It is 
found that the first (hydrogen-bonded) isomer has a slightly smaller dipole moment than the 
pure MgSOq molecule (12.15 Db vs. 12.85 Db), while the second isomer has a much 
larger dipole moment (17.30 Debye). Binding energies of the two MgSC>4 • H2O 
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complexes, calculated at the MP2/TZP level of ab initio theory, are 42 and 41 kcal/mol. All 
energetics are summarized in Table 1. 

Two isomers are also found for the complex of MgSC >4 with two water molecules 
(Fig.3). The lowest energy isomer has one water molecule that is connected to the Mg atom 
only, and the other water molecule that in addition forms a hydrogen bond with an oxygen 
atom from MgSC> 4 , similar to the hydrogen bond in the complex with one water molecule 
(Fig. 2a). This isomer has a dipole moment of 13.32 Debye. Its total binding energy is 79 
kcal/mol or 39.5 kcal/mol per water molecule. The second isomer has a “bridge” of two 
water molecules connecting the Mg atom with one of the oxygens of the SO 4 group. This 
isomer is about 12 kcal/mol higher in energy than the first (lowest energy) one. The two 
hydrogen bonds of the “bridge” in the second isomer are also very strong, but the 
participating O-H bond lengths are not as elongated as in the 3a and 2a structures. 

The two isomers of the complex of MgS04 with three water molecules are shown in 
Fig. 4. Again, the lowest energy isomer has only one hydrogen-bonded water and two “free” 
water molecules bound only to Mg atom. The total binding energy of this complex is 109 
kcal/mol or 36.2 kcal/mol per water molecule, and its dipole moment is 13.27 Debye. The 
hydrogen bond in this complex is similar to hydrogen bonds in the lowest energy complexes 
of MgS04 with one and two water molecules, but is slightly weaker. This hydrogen bond is 
the strongest in the 1 : 1 complex and becomes somewhat weaker as the number of water 
molecules increases. The second isomer has a two-water "bridge" and a hydrogen-bonded 
water molecule. It is 4.6 kcal/mol higher in energy than the other (most stable) isomer and has 
a smaller dipole moment of 1 1.30 Debye. The hydrogen bonds of this isomer are not as 
strong as in other isomers, including the hydrogen bonds of the “bridge”, which are weaker 
than those in the similar bridge of the complex with two water molecules (structure 3b). 
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B. Vibrational spectra 

Harmonic and anharmonic (CC-VSCF) vibrational frequencies as well as infrared and 
Raman intensities for the isolated water molecule, isolated MgS04, and the complex 

MgS0 4 1K,0 obtained in this study are presented in Table 2. Comparison of vibrational 

frequencies of the water molecule in the 1:1 complex with MgSOa (where it forms the 
hydrogen bond with one of the oxygen atoms of SO4) with vibrational frequencies of the pure 
water molecule, shows that the frequency of one of the O-H stretching vibrations is reduced 
by about 1400 cnr 1 at the harmonic level (down to 2450 cm' 1 ). When anharmonic corrections 
are included, this frequency shift due to complexation of H 2 0 with MgS0 4 is even more 

significant, almost 2000 cm'1 The CC-VSCF frequency for this O-H stretch is 1642 cm - ', 
which is in the region of water bending frequencies. However, its infrared and Raman 
intensities are very high, much higher than the corresponding intensities of the water bending 
vibrations. Therefore, this unusually red-shifted O-H stretching vibration should be 
observable and easily distinguishable from bending vibrations. 

Since the anharmonic correction for this O-H stretch is unusually large (about 800 cnr 
!), we have checked the stability of this result by computing more dense grids for the pair 
coupling potentials of this mode with all other normal modes of the complex. 16x16 point 
grids were computed and interpolated into more dense 32x32 point grids. This changed the 
vibrational frequencies of the complex by only a few cm’ 1 , including the value of the red- 
shifted O-H stretch that changed only slightly, from 1642 to 1639 cnr 1 . This shows that the 
calculated frequency for this unusually red-shifted and anharmonic vibration is not a numerical 
artifact, but a stable result. 

It can be also seen from Table 2 that one of the S=0 stretching vibrations of MgS04 is 

also red-shifted, but to a lesser extent, by about 200 cm 4 upon complexation with one water 
molecule, and its infrared intensity is significantly enhanced. In the second isomer of 
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MgS0 4 H 2 0 (Fig. 2b), the water molecule does not form a hydrogen bond, and its 

frequencies are not changed as much as in the hydrogen-bonded isomer relative to the isolated 
water molecule. However, its infrared intensities are also significantly enhanced and therefore 

should be distinguishable from the pure water vibrations. 

The effects found for the complex of MgS0 4 with one water molecule are also present 
in complexes of MgS0 4 with two and three water molecules. The O-H stretching frequencies 
of these complexes are significantly shifted to the red compared with the isolated water, but to 
a somewhat lesser extent than in the 1 : 1 complex. Note that anharmonic corrections were not 
obtained for the complexes of MgS0 4 with two and three molecules, since CC-VSCF 
calculations for such large systems are very expensive. However, comparison of harmonic 

frequencies for these systems with those of the 1 : 1 MgS0 4 -H 2 0 complex shows that the O-H 
stretching frequencies of the larger complexes are red-shifted to a lesser extent. The red 
shifted O-H stretches of the MgS0 4 2H z O complexes are located at about 3400, 2800, and 

2500 cm' 1 at the harmonic level, while those of the MgS0 4 -3H 2 0 complexes are at about 

3450, 2940, 2800, and 2660 cm' 1 (harmonic level). We expect these frequencies to go down 
significantly when anharmonic corrections are included. In our estimate, the anharmonic 
corrections should be of the order of several hundreds cm 1 , 350 to 600 cm on average. 
Thus, the red shifted O-H stretches should be observed in the region between 3200 and 2000 
cm' 1 . This region is very different from pure water and water cluster frequency regions: 3600- 

3800 cm' 1 for O-H stretches and 1600 cm' 1 for HOH bends. The infrared intensities of these 
vibrations are about 20-30 times higher than those of the strongest pure water O-H vibrations. 
Observation of such intense infrared vibrations in the 2000-3200 cm' 1 region should be 


indicative of the presence of hydrated salts. 
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In addition to O-H stretches, water bending vibrations also change upon complexation 
with magnesium sulfate. Their frequencies are somewhat higher (blue shifted) in the 
complexes (up to 1660-1710 cm' 1 ) relative to the pure water molecule (about 1600 cm ), and 

they also have larger infrared intensities (see Tables 2-4). 

Unlike O-H stretches for water molecules, the vibrational frequencies of the MgS0 4 
fragment are found to be nearly harmonic, except for the S-0 stretching vibration of the S=0 
bond participating in the hydrogen bonding with water. This S=0 stretching vibration goes 

down from 1009 to 935 cm' 1 in the MgS0 4 H 2 0 complex when anharmonic corrections are 

included. The highest S=0 stretching frequency is located around 1310-1320 cm , and its 
position almost does not change in the different complexes. This vibration, and the other S=0 
stretch in the region of 1000 cm' 1 , have very large intensities and can serve as unique 
fingerprints of the presence of sulfate fragments. 

Conclusions 

Geometrical structure, properties, and vibrational spectra of pure and hydrated 
MgS04 salts are calculated using electronic structure theory at the second order Mpller- 
Plesset perturbation theory (MP2) level with a triple-zeta + polarization (TZP) basis set. For 
smaller systems (pure MgS04 and its complex with one water molecule) anharmonic 
corrections to vibrational frequencies are computed using the correlation corrected vibrational 
self-consistent field (CC-VSCF) method. 

Very significant differences are found between vibrational spectra of water molecules 
in complexes with MgS0 4 and pure water. Some of the water O-H stretching frequencies are 
shifted to the red very significantly upon complexation with magnesium sulfate, and their 
infrared intensities are greatly enhanced. This effect is due to very strong hydrogen bonds 
formed between water molecules and the oxygen atoms of MgS04- The red shift is largest 
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for the smallest 1:1 MgS0 4 H 2 0 complex, where the O-H stretching vibration is shifted by 

almost 2000 cm’ 1 and is located at about 1640 cm' 1 . This vibration should be easily 

distinguishable from HOH bending vibrations since it is calculated to have a very high 
infrared intensity (almost 20 times that of the HOH bending in the pure water). Strongly red 
shifted O-H stretches are also present in complexes of MgS0 4 with two and three water 
molecules, but this effect becomes less significant as the number of water molecules around 
MgS04 increases. However, even for the situation where MgS04 has high number of 
hydration waters, the hydrogen bonded O-H stretching frequency should be observed in the 
region between 2000 and 3200 cm’ 1 , which differs significantly from pure water and water 

cluster frequency regions (3600-3800 cm' 1 ). Water bending frequencies in complexes may 
be easier to identify, since they do not change as much as the stretching frequencies. They are 
found to shift to the blue by up to 100 cm’ 1 upon complexation with magnesium sulfate. In 

addition, the S0 2 stretching vibrations are found at lower frequency regions (900-1300 cm l ) 

than the water vibrations. They can serve as unique identifiers for the presence of sulfate 
salts. The predicted infrared and Raman spectra should be of valuable help in the design and 
analysis of experimental observations of the icy surface of Jupiter’s moon Europa that 
possibly contains hydrated MgS0 4 salts. 

For the 1:1 MgS0 4 H 2 0 complex, the anharmonic correction of the red-shifted O- 

H stretch, calculated using CC-VSCF, is very large, about 800 cm’ 1 . Anharmonic 

corrections of this size have not been obtained previously for neutral hydrogen-bonded 
complexes. Experimental measurement of this transition will be very useful as a test of our 


theoretical treatment. 
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Table 1. Total energies (hartrees) and dipole moments (Debye) 
of complexes of MgS04 with water. 


Total energy 


MgS0 4 

H2O ~~ 

MgS0 4 lH 2 0 isomer a 
isomer b 
MgS0 4 -2H 2 0 isomer a 


-897.611838 

-76.274033 

-973.953258 

-973.951329 

-1050.285492 


Dipole moment 
12.846 
1232 
1L66T 

17.308 

13.318 


isomer b 


MgS0 4 -3H 2 0 isomer a 


-1050.266681 

-1126.607081 


13.249 

13.272 


isomer b 


-1126.599753 


11.304 
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Table 2. Vibrational frequencies (cm 1 ), infrared (Debye 2 / A 2 amu) and 
Raman (A 4 / amu) intensities of MgSC> 4 , H 2 0, and MgS0 4 - 1H 2 0 complex. 



frequencies 

intensities 


harmonic 

CC-VSCF 

IR 

Raman 

h 2 o 





asym. OH stretch 

4015 

3797 

1.19 

34.99 

sym. OH stretch 

3877 

3691 

0.17 

72.48 

HOH bend 

1603 

1548 

1.80 

6.21 

MgS 04 





asym. S=0 stretch 

1347 

1326 

7.01 

22.13 

sym. S=0 stretch 

1167 

1153 

6.33 

15.67 

MgS0 4 lH 2 0, isomer a 





free OH stretch 

3958 

3802 

3.16 

83.49 

bonded OH stretch 

2451 

1642 

33.97 

51.36 

HOH bend 

1676 

1547 

7.07 

7.39 

free S=0 stretch 

1333 

1318 

8.90 

11.44 

bonded S=0 stretch 

1009 

935 

10.95 

13.06 

MgS0 4 lH 2 0, isomer b 


| 



asym. OH stretch 

3946 


5.41 

26.45 

sym. OH stretch 

3845 


3.17 

64.97 

HOH bend 

1685 


3.36 

1.03 

asym. S=0 stretch 

1330 


7.61 

14.50 

sym. S=0 stretch 

1162 


6.98 

14.32 
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Table 3. Vibrational frequencies (cm- 1 ), infrared (Debye 2 / A 2 amu) and Raman (A 4 / amu) 
intensities of complexes of MgSC >4 with two water molecules 



harmonic 

frequencies 

Intensities 
IR Raman 

MgS0 4 -2H 2 0, isomer a 




W 1 : free OH stretch 

3961 

2.59 

81.04 

bonded OH stretch 

2499 

35.30 

63.47 

HOH bend 

1693 

5.22 

8.14 

W2: asym. OH stretch 

3964 

4.78 

51.83 

sym. OH stretch 

3808 

2.45 

54.84 

HOH bend 

1655 

4.32 

1.39 

free S=0 stretch 

1320 

9.00 

9.04 

bonded S=0 stretch 

1020 

10.80 

12.62 

MgS0 4 -2H 2 0, isomer b 




W 1 : free OH stretch 

3973 

4.19 

83.27 

bonded OH stretch 

2807 

33.38 

72.44 

HOH bend 

1705 

4.95 

4.65 

W2: asym. OH stretch 

3949 

2.83 

60.03 

sym. OH stretch 

3419 

21.69 

125.85 

HOH bend 

1667 

1.61 

4.17 

free S=0 stretch 

1337 

7.74 

11.25 

bonded S=0 stretch 

1163 

6.85 

13.35 
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Table 4. Vibrational frequencies (cm 4 ), infrared (Debye 2 / A 2 amu) and Raman (A 4 / amu) 
intensities of complexes of MgSOa with three water molecules 


harmonic 

frequencies 


Intensities 
IR Raman 


MgS0 4 -3H 2 0, isomer rd 


W 1 : free OH stretch 


3947 


2.17 


69.88 


bonded OH stretch 


2664 


33.75 


75.78 


HOH bend 


1710 


2.97 


5.05 


W2: free OH stretch 


3954 


3.74 


45.74 


bonded OH stretch 


3770 


2.62 


68.30 


HOH bend 


1646 


2.89 


3.31 


W3: free OH stretch 


3966 


4.11 


55.33 


bonded OH stretch 


3757 


2.63 


51.67 


HOH bend 


1628 


4.10 


3.12 


free S=0 stretch 


1312 


8.58 


8.01 


bonded S=0 stretch 


1040 


10.26 


11.39 


MgS0 4 -3H 2 0, isomer b 


W 1 : free OH stretch 


3958 


2.66 


82.81 


bonded OH stretch 


2806 


24.71 


63.79 


HOH bend 


1699 


4.44 


0.75 


W2: free OH stretch 


3971 


3.73 


76.82 


bonded OH stretch 


2941 


35.38 


94.23 


HOH bend 


1711 


2.90 


5.08 


W3: free OH stretch 


3951 


2.90 


53.09 


bonded OH stretch 


3449 


21.47 


129.28 


HOH bend 


1662 


2.38 


3.77 


free S=0 stretch 


1317 


8.48 


8.26 


bonded S=0 stretch 


1067 


9.00 


7.65 
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Figure 1. Equilibrium structures (angs.) and natural charges of pure MgS0 4 and water molecules. 
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Figure 2. Equilibrium structures (angs.) and natural charges of complexes of MgS0 4 with one water 
molecule. Structure (a) corresponds to the lowest energy isomer. 





Figure 3. Equilibrium structures (angs.) and natural charges of complexes of MgS0 4 with two water 
molecules. Structure (a) corresponds to the lowest energy isomer. 
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Figure 4. Equilibrium structures (angs.) and natural charges of complexes of MgS0 4 with three water 
molecules. Structure (a) corresponds to the lowest energy isomer. 





